arXiv:1508.06348vl [cond-mat.supr-con] 26 Aug 2015 


Resonance in the nonadiabatic quantum pumping of the 
time-dependent Josephson junction 


Rui Zhu* and Mi Liu 

Department of Physics, South China University of Technology, 

Guangzhou 5106fl, People’s Republic of China 

Abstract 

In this work, we investigated the nonadiabatic transport properties of the one-dimensional time- 
dependent superconductor-normal metal-superconductor (SNS) Josephson junction biased by a 
current source and driven by a high-frequency-ac-gate-potential applied to the normal-metal layer. 
BCS superconductors are considered and treated with the time-dependent Bogoliubov-de Gennes 
equation. Using Floquet theory, we compute the transmission coefficients and the Wigner-Smith 
delay times as a function of the incident energy and find that they display resonances when one 
of the electron or hole Floquet wavevectors coincides with the bound quasiparticle state within 
the superconducting energy gap. The resonance varies with the phase difference between the two 
superconductors as a result of the bound quasiparticle level displacement. The supercurrent flowing 
through the SNS junction is dramatically enhanced by the resonances. 
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I. INTRODUCTION 


The scattering process often involves interferences between different quantum paths, in 
which constructive interference corresponds to resonance and destructive interference to 
antiresonance of the transmission. Asymmetric antiresonance with a minimum followed 
by a maximum is called a Fano resonancei*^. There have been a great number of studies 
devoted to resonance and/or antiresonance in various quantum processes, such as scattered 
by an Anderson impurity^, tunneling through a quantum dot^'-, scattering from a donor 
impurity in an electron waveguide^*^, transport in spin inversion devices^, Mie scattering in 
plasmonic nanoparticles and metamaterials^, and etc. In nonadiabatic quantum pumping, 
Floquet sidebands are formed by high-frequency oscillating potentials. In the case of a time- 
dependent quantum well, Fano resonance occurs in the transmission spectrum^^— , pumped 
shot noisei^ii^, and Wigner-Smith delay times (WSDT/IT'I^, when the energy or wavevector 
of one of the Floquet levels matches the quasibound level inside. To our knowledge, these 
resonant effects have not been discovered in the time-dependent superconductor-normal 
metal-superconductor (SNS) Josephson junction driven by a nonadiabatic electric potential 
applied to the normal region. 

By means of the Josephson effect, supercurrents can flow through the SNS junction biased 
by a current source!^. Bound quasiparticle states exist in the normal region of the SNS junc¬ 
tion, each consisting of equal probabilities of particle and hole states. The Josephson effect is 
related to the Andreev reflection, the latter of which is defined by reflection of a particle into 
a hole at a pair-potential boundary with no change of current^. It can be conjectured from 
previous studies that in nonadiabatic quantum pumping of the SNS junction, resonance is 
possible to occur in the supercurrent when one of the Floquet sidebands coincides with one 
of the bound quasiparticle states within the energy gap of the superconductor, which will 
be conhrmed by the present theoretic work. 

Our approach is based on various development in methodology and concept recently. In 
classical mechanics, the duration time of a scattering process can be defined by the energy 
derivation of the action. In the absence of a time operator in quantum mechanics, the 
quantum analog of the duration time can be defined as the energy derivative of the quantum 
mechanical phase shift r = hd(f)/dE during the scattering process. In a multi-channel or 
dynamic scattering process, the distribution of the WSDT ri, • • • ,tn are the eigenvalues of 
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the N X N Wigner-Smith matrix Q{E) = —ihS^dS/dE with S the scattering or Floquet 
scattering matrix— and N = 4 (2?7,jnax + 1) with the maximal Floquet channel index 
in the present dynamic electron-hole system. The density of states (DOS) p {E) is directly 
related to Q{E) and the WSDT by 

p (E) = (27rh)-'Trg (E) = {2nh)-^ (1) 

n 

The parametric conductance derivatives of a quantum dot was investigated by relating it to 
the distribution of the Wigner-Smith time-delay Matrix—. And also the WSDT were inves¬ 
tigated in the scattering by strong time-periodic driving helds^'^. Experiments can yield 
information on the distribution of the WSDT. The driving held can be realized by applying 
an ultrahigh intensity laser or a local ac-signal top gate. In the research of the charge turn¬ 
stile based on the superconducting hybrid structure, it was found that quantization of the 
current is affected by Andreev rehection and Cooper pair-electron cotunnelingi^. Nonadia- 
batic effects in braidings of Majorana fermions in topological superconductors were studied 
using the time-dependent Bogoliubov-de Gennes (BdG) equations^S. The time-dependent 
BdG equation can be used to describe time-evolution of BCS superconductors under a 
parametrically time-dependent Hamiltoniaii22>^. In this paper, we consider the transport 
properties of the time-dependent SNS Josephson junction biased by a current source and 
driven by a high-frequency-ac-gate-potential applied to the normal-metal layer. We use the 
Floquet scattering matrix method to solve the time-dependent BdG equation and relate the 
Floquet scattering matrix with the WSDT and the DOS. With the general relation between 
the Josephson supercurrent and the quasiparticle excitation spectrunt^S as well as the bound 
quasiparticle states and the DOS of the continuous spectrum, the supercurrent is calculated. 


II. THEORETIC FORMALISM 

We consider the nonadiabatic pumping properties in the one-dimensional time-dependent 
SNS Josephson junction biased by a current source and driven by a high-frequency-ac-gate- 
potential applied to the normal region. The time-dependent electric potential has the form 
of Vicosut. Width of the normal region is L. The considered model is sketched in Fig. 
1. The driving held can be realized by applying a local ac-signal top gate. Assuming the 
junction is located in the x-direction, the time-dependent BdG equation can be expressed 
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The time-dependent BdG Hamiltonian is 
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( 2 ) 


with 


-f^BdG {t) — 


h {x, t) = 


h{x,t) A{x) 
(a;) —h{x,t) 


(3) 




2,ndx‘ + 

(4) 

(a;) 9 {L — x) Vi cos (cat), 

(5) 

f 

Aoe'^/^, X < 0, 


0 , 0 < X < L, 

(6) 

Aoe“'^/^, X > L. 



A (a;) = < 


9{x) is the step function. 

In advance of the time-dependent treatment, we consider the bound quasiparticle states in 
the normal region within the energy gap of the superconductor. The quasiparticle energy e: is 
measured with respect to the Fermi energy Ep- The current-flux normalized eigenfunctions 
of the static BdG equations are^^ 


with 


^N,e = 


'4’nx = 




rei’si,e + rhi^si,h^ a: < 0, 

+ HN,e + C^N,h + 0 < X < L, 

tei^S2,e + th‘lps2,h^ X > L, 

(fee) exp {±ikeX) , 

(4)”^^%xp {±ikhx ), 

{2qe,h)~^^‘^{£‘^/Al - l)"^'^^exp {±iqe,hx) , 


(7) 


( 8 ) 



while 4’s2eh V^e,h replaced by ri2e,h and exp {±iqe,hx) replaced by exp [±ige,/i {x — L)]. 

r 1 1/2 

Here, = {2m/h^f^‘^{EF + qe,h = {2m/h^f^'^ Ep (Te,h{e^ - , r]^e,h = 

0/2-|-(Te,ft arccos (e/Ao), r]2e^h = —0/2-|-cTe,h arccos ( 5 /A 0 ), <Je = 1, and ah = —1. The square 
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roots are to be taken such that Rege,/i > 0, Imge > 0, and Ivnqh < 0. The function arccost G 


for t > 1. Solvability of the 


[0,7r/2] for 0 < f < 1, while arccost = —iln t+ — 1)^^^ 
continuity equations of the wave functions and their derivatives at a; = 0 and x = L gives 
rise to the secular equation. Roots of £ of the secular equation are the bound quasiparticle 
energies measured from the Fermi energy Ep. These bound states consist of equal 
probabilities of particle and hole states. Numerical results of Ef, as functions of Aq and 0 
are shown in Fig. 2. For a small superconductor energy gap Aq relative to the Fermi energy 
Ep and a small normal-region width L considered here, only a single bound quasiparticle 
state is sustained in the normal region^. The bound quasiparticle level is close to the gap 
surface and increases linearly with the gap width. It varies in an absolute trigonal function 
with the phase difference 0 between the two superconductors. It approaches maximum at 
0 = 0 and 271 and approaches minimum at 0 = tt. These results reproduce those of Bardeen 
and Johnsoii^ in 1972. 

Now we use the Floquet scattering theory to solve the time-dependent BdG equation ([2]). 
Wave functions in the superconductor and normal regions can be written as 


+ 00 




^ — i£nt 


0si,n, a; < 0, 

0 < a: < L, 

0S2,n, X> L, 


(9) 


with En = s + nhu, e an energy within the static continuous quasiparticle spectrum above 
the energy gap, n an integer varying from —oo to -|-oo, and 


051,n = 0^051,e,n + ^L0Sl,ft,n + Cn051,e,n + 4051, 

+ OD 

0V,n = X] [(®m0V,e,m ^m'4’N,e,ra “I" '^m0V,h,m) '^n-m (hl/hoi)] , (10) 

m=—oo 

052,n = <002,e,n + <0S2,/^,n + 4052,e,n + 4002,h,n- 

Here, a4 and bn'~ are the probability amplitudes of the incoming electron and hole waves 
from the left/right superconducting electrode of the nth Floquet channel with energy Sn, 
respectively, while c4 and dn^ are those of the outgoing electron and hole waves. And we 
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have 


^N,e,n 
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( 11 ) 


(2ge,h;n) - 1) "^^exp {±iq^^h.nX) , 


-1/4, 


while e ft-n 1^^® V^e,h;n replaced by i/2e,h;n and exp (±ige,/i;na^) replaced 

by exp [±iqe^h-,n {x - L)] with ke,h-,n = {2m/h‘^Y^‘^{EF + ae,h£nY^‘^, qe,h;n = 

{2m/h‘^f^‘^ Ep + (Je,h{£l- , hle,h;n = 4^/2 + ae,ft arccos (£„/Aq), and 
h2e,/i;n = “0/2 + (Je,h arccos (e„/Ao). Jn{x) are the nth-order hrst kind Bessel func¬ 
tions. 

By continuity of 0 and dip/dx at the normal-superconducting interfaces, the electron-hole 
Floquet scattering matrix expressed as 


^ Cn ^ 



= E 


m 

< 

m 

\<) 


[bln) 



/ r, 
r, 
t 

\k 


ee 

j^eh 

.fee 
^ nm 

,/eh 
^ nm 

\ 

nm 

nm 

he 

r^hh 
' nm 

,fhe 
^ nm 

,fhh 
^ nm 


nm 


ee 

j-eh 

^nm 

tee 

feh 


nm 

^ nm 

^ nm 


he 

fhh 

^nm 


fhh 
^ nm 

/ 

nm 

' nm 


/ aL \ 

6: 


fj. 

-^m 


( 12 ) 




can be obtained by matrix algebra (see the Appendix). We dehne the total electron and 
hole reflection and transmission coefficients as 
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(13) 


The summation is over all propagating modes and with a cutoff of nmax in numerical treat¬ 
ment. The Wigner-Smith matrix can be obtained by Q (e) = —ihS^dS/de. The DOS p{e) 
is directly related to Q (e) and the WSDT by Eq. ([T]). We can also dehne the total 
WSDT as Tws = Biased by a current source howing from the Si to the S 2 super- 

n 

conducting electrodes and driven by a high-frequency-ac-potential in the normal region, the 
supercurrent howing from the Si to the S 2 region can be calculated by the general relation^^ 


I = h + /2, 


(14) 
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with 



(15) 


(16) 


III. NUMERICAL RESULTS AND INTERPRETATIONS 

In our numerical treatment, experimentally realistic parameters are used. Ep = 20 meV, 
Ao = 2 to 3 meV, ho; = 5 to 6 meV, Vi = 1 meV, and Umax = 1 > Vi/IEj. Numerical 


results of the total electron and hole reflection and transmission coefficients dehned in Eq. 


f[T3|l are shown in Fig. 3. Sharp resonances occur in Ree, Reh, and Teh', sharp antiresonances 
occur in Tee, when one of the electron or hole Floquet sidebands coincides with the bound 
quasiparticle state Ef, by Ee — hu = E^ in the electron channel and hw — Eh = Ef, in the 
hole channel within an accuracy of 0.001 meV. Eh = 1.99519 and 1.99487 meV can be 
obtained by solution of the secular equation for 0 = 0 and 0.025 in radian, respectively. 
Current conservation secures unitarity of the Floquet scattering matrix S. It can be seen 
from Fig. 3 that Ree + Reh + Tee + Teh = 1, which conhrms unitarity of the S matrix. 
From Fig. 2 it can be seen that Eh decreases as 0 in creases for 0 < tt. As 0 increases 
from 0 to 0.025 in radian, resonance from the electron Floquet channel occurs at a larger 
energy and resonance from the hole channel occurs at a smaller energy according to the 
relation Ee — hw = Eh and hw — Eh = Eh- It can also be seen that the resonance peaks of 
electrons are broader than that of the holes. This is because that the virtual mass of the 
hole is larger than that of the electron and the characteristic time of the former is longer 
than the latter, which is prominently seen in Fig. 4. In the range Aq < e < ho; + Aq we 
considered, ke^h-,o, ke,h-,±i, QeAfi, and ge,h;±i are all real in the normal region; ke^hfi, keA,+i, 
QeAO, and q'e,/i;+i are real in the superconducting regions; ke^-i is imaginary and QeA-i 
complex in the two superconducting regions evanescent to ±cxd. The incident quasiparticle 
state transfers through the propagating quasiparticle Floquet sideband and interferes with 
the direct transmission giving rise to a resonance. 

Sharp resonances also occur in the r-ws (see Fig. 4) featuring the characteristics of 
the Floquet scattering matrix. The resonant rws of the hole is larger than that of the hole 
because of larger virtual hole mass, which is natural in electron-hole systems. Our numerical 
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results also show that r-ws outside of the resonant peak is to 10“^^ s, a normal WSDT 
in nanoscale tunneling. Values of the rws at resonance are as high as 10“^ s, which is a 
dramatic enhancement in comparison with the normal 10“^^ s and possible to be observed 
in experiment. 

The energy gap of BCS superconductors Aq is approximately related to the temperature 
by Aq = 3 . 2 fc 5 Tc(l — Twith Tc the critical temperature^^. From Fig. 2 it can be seen 
that the bound quasiparticle energy increases linearly with Aq for hxed Ep, L, and (j). It can 
also be seen from Figs. 3 and 4 that resonances in the reflection and transmission coefficients 
and Tws have extremely sharp peaks with their widths smaller than 0.0002 meV in a wide 
energy window of 5 meV. The peak energy is Eg — hw = Eb and hw — Eh = Eb within an 
accuracy of 0.001 meV. These results based on the time-dependent SNS Josephson junction 
suggest a potential way to detect the bound quasiparticle state energy, the superconducting 
energy gap, and hence the superconducting transition temperature. 

Supercurrent of the considered device are calculated by Eqs. cai to fllOp . The numerical 
results are given in Fig. 5. The supercurrent flowing through a static Josephson SNS 
junction is^^ in the order of meV • e/h. It can be seen from Fig. 5 that the supercurrent 
is dramatically enhanced in magnitude by the nonadiabatic driving force. This can be 
interpreted by the change in the continuous DOS above Aq as a result of the Floquet states 
formed. As demonstrated in Eq. ([T]), the DOS is proportional to the sum of the eigenvalues 
of the Wigner-Smith matrix Q (e) = —ihS^dS/dE with S the Floquet scattering matrix. The 
eigenvalues of the Wigner-Smith matrix corresponds to the WSDT. The Floquet sidebands 
are formed in the high-frequency driven transport processes. The sidebands contribute to 
the dynamic DOS and enhance the supercurrent. From Figs. 3 and 4, it can be seen 
that when one of the Floquet sidebands coincides with the bound quasiparticle state, sharp 
resonances occur in the transmission and WSDT spectrum, which further enhances the 
supercurrent. Enhancement by the resonances is vital as the direction of the supercurrent 
is reversed against the original current bias. Energies of the peaks in rws and hence in p 
corresponding to the hole-Floquet-channel resonance decreases with 0 while energies of the 
peaks corresponding to the electron-Floquet-channel resonance increases with 0 governed 
by the relations Eg — hw = Eb and hw — Eh = Eb as Eb decreases with 0 for 0 < tt. Ii in 
Eq. flT^ is proportional to —dEb/dcf) and is positive for 0 < vr and negative for 0 > vr. I 2 in 
Eq. ffTHj) is proportional to the integral of —dp/dcf) over e. Contribution of the hole-channel 
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peak in p to 1-2 is negative for 0 < tt and positive for 0 > tt, while contribution of the 
electron-channel peak in p to I 2 is reversed. Since the hole-channel peak in p is higher and 
sharper than the electron-channel, contribution of the former to I 2 overweighs that of the 
latter, giving rise to the sign reversal of the supercurrent. As the resonances are extremely 
strong, the supercurrent in the nonadiabatic process is several orders larger than the static 
supercurrent. Since —dp/d(j) is nonzero even when p is small at the energies away from the 
resonant peaks, the supercurrent varies in an irregular pattern as it is an integral result 
of —dp/d(j). For (p > n, increases with 0. The supercurrent is positive or negatively 
small as a combined result of the resonant peaks and the integral of —dp/dp over the large 
range of energy. The resonant peaks are sharper for stronger driving forces. As a result the 
super current for hw = 6 meV is larger than hw = 5 meV for p < 71 . 


IV. CONCLUSIONS 

In conclusion, resonances are observed in the nonadiabatic transmission and WSDT spec¬ 
tra of the time-dependent Josephson SNS junction as a result of quantum path interferences. 
In the picture of Floquet scattering, the electron and hole Floquet sidebands are formed. 
When one of the electron or hole Floquet sidebands coincides with the bound quasiparticle 
state within the energy gap of the superconductors at Eg — hoj = Ei, and tvjj — Eh = Eb, 
resonances occur in i?ee, Reh, Teh, and rws; antiresonances occur in Tee- The resonances are 
extremely sharp suggesting a potential experimental determination of Eb- Using the Wigner- 
Smith matrix, the dynamic DOS and supercurrent biased by a current source are obtained. 
Sharp resonances in the dynamic DOS consisting of Floquet modes greatly enhance the 
supercurrent in comparison with the condition without a pumping force. 
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VI. APPENDIX: DERIVATION OF THE FLOQUET SCATTERING MATRIX 


Continuity equations of xjj and d'lfj/dx [The Floquet state wave functions are defined in 
Eqs. ([9]) to ffTTjl ] at the normal-superconducting interfaces a: = 0 and x = L can be expressed 
in matrix form as follows. 

MtJ + + MlJ + Mid' = M‘,a + M‘,c, (17) 

Mi>‘ + M’lJi' + + M^jd' = Mj6 + Mjd, (18) 

A4%o + Ml_c = + Mlb' + Aqy + Miy, (19) 

Mlb + M^d = MJ.o' + + M^d--, (20) 

+ Ml,J + Ml,_c‘ + = A4%a + Af^c, (21) 

+ M‘u^J + A4,+d' = Mlb + Mid, (22) 

Mia + Mlc = Ml,_a'' + Ml,^V + Af.'^^^.c’' + Ml^_dd, (23) 

Mlb + Mid = + Ml.^b’- + A/L.^c’- + A/ij.tT. (24) 

a^/’’, 1)*'^’’, and db^’’ are column vectors made up of elements a[l, bl, cl, and dl, 
respectively. Elements of the coefficient matrices are dehned as follows: 

= exp (2r/le,ft.„/2) (25) 

\ / nm 

(Ml'f) = exp {-ir]le,h-n/‘^) (2ge,/i;n)”^'^^(4/2^0 - (26) 

V / nm 

(Mf) = Jn-m {Vjhw) (27) 

(M 2 u) = exp {ir]2e,h-,n/2) (2ge,h;n)"^^^(4/^0 - (28) 

) = exp (-* 772 e./i;n/ 2 ) {2qy^){el!- l)(29) 

(m^^) = exp{±ike,hirnL) {ke,h-,m)~^^‘^Jn-m{Vi/hw) , (30) 

= ±iqe,h-,n exp (^77le,/^;r^/2) (2ge,ft.„)(4/2^0 - l) (31) 

(^4dl±) = ±iqe,h-,nexp {-irile,h-,n/y (2ge,/i;n)"^^^(4/^0 “ (32) 

(mI'£) = ±ike^h-,m{ke,h-,m)~^^‘^ Jn-m {Vi/ku) , (33) 

(m1u2±) = exp {ir]2e,h-,n/y (2ge,/i;n)"^^^ (4 /^^0 “ l) (34) 

V / nm 
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m: 


e,h 


4 d 2 ± 


= ±ige,h;nexp (-i772e,ft.„/2) (2ge,ft.„) - l) 

^S± ) ~ i^^e,/i;m6Xp (A;e,/i;m) ^ Jn—m • 

/ nm 

The matrix equations ffTTl) to fl2T|) can be transformed into larger matrix equations 

/a'\ ( Mj 0 Mj 0\/a\ 


nm 
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( 36 ) 
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Then we have 


^ a'’ ^ 


/a'\ 

6" 

= (Mr)“^ • (M02) ■ (MOl)"^ • {Ml) • 
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d 



\d) 


Uv 


By dehning 


Mltor — 


Mitor = {Mr)~^ ■ (M02) • (M01)“^ • (Ml ), 

/ W^a Wab Wac \ 

Wba Wbb Wbc Wbd 
W,a W,b PBcc W,d 

\ Wda Wdb Wdc Wdd J 

( -W^c -Wad 0 0 \ 

—Wbc —Wbd 0 0 
-Wcc -Wed 1 0 

V —Wdc —Wdd 0 1/ 


Med = 


(37) 
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(43) 


Mnh = 


/ Waa -1 0 \ 

Wha Wbb 0 —1 

Wca Web 0 0 

V Wda Wdb 0 0 y 

we can obtain the electron-hole Floqnet scattering matrix as 


S = (Mrf)-' . 


(44) 


In Eqs. fH2]) and fH5]l . “0” is a zero (2nmax + 1) x (2nmax + 1) rnatrix and “1” is a nnitary 
(2nmax + 1) X (2ninax + 1) niatrix with rimax the maximal Floqnet channel index. 
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FIG. 1; Model for the one-dimensional time-dependent Josephson SNS junction. Width of the 
normal region is L. Order parameters of Si and S2 regions have the same absolute value of Aq and 
phase difference of (j). When Aq and the Fermi energy Ep is small, only one bound quasiparticle 
state with the energy of Ei, is confined in the normal region consisting of equal probabilities of 
particle and hole states. A high-frequency-ac-gate potential of Vi cos cut is applied to the normal 
region. Driven by the time-dependent potential, Floquet sidebands are formed with energy spacing 
hui. When the continuous spectrum above the energy gap has the electron or hole Floquet sideband 
falling within the energy gap and coincident with the bound state, resonance or antiresonance occurs 
in the transmission coefficients and the Wigner-Smith delay times. 
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FIG. 2: Bound quasiparticle energy Eh measured from the Fermi energy Ep. Its variation as func¬ 
tions of Aq and 4> corresponds to black and red symbols, respectively. Corresponding coordinates 
are in the same colors as the symbols. The parameters are Ep = 20 meV, L = 10 A, iji = 0 for 
variation as a function of Aq and Ao=2 meV for variation as a function of </>, respectively. 
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FIG. 3; Numerical results of the total reflection and transmission coefficients (a) Ree, (b) Reh, 
(c) Tee, and (d) Teh as a function of the incident quasiparticle energy e for different 4>. Sharp 
resonances occur in Ree, Reh, and Teh and sharp antiresonances occur in Tee, when one of the 
electron or hole Floquet sidebands coincides with the bound quasiparticle state Eh hy Ee — fiu} = Eh 
in the electron channel and hu — Eh = Eh in the hole channel within an accuracy of 0.001 meV. 
Insets are zoom-in of the two resonant peaks. The parameters are Ep = 20 meV, T = 10 A, Aq = 2 
meV, hu} = 5 meV, Vi = 1 meV, and Umax = 1- Eh = 1.99519 and 1.99487 meV for (;/i = 0 and 
0.025 in radian, respectively. 
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FIG. 4: Numerical results of the total WSDT rws as a function of the incident quasiparticle energy 
e for different cj)- Sharp resonances occur when one of the electron or hole Floquet sidebands 
coincides with the bound quasiparticle state E), hy eg — fiuj = E), in the electron channel and 
huj — Eh = Ef, in the hole channel within an accuracy of 0.001 meV. Insets are zoom-in of the two 
resonant peaks. The parameters are the same as Fig. 3. 
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FIG. 5: Supercurrent as a function of cj) for different driving frequencies. The parameters are 
Ep = 20 meV, L = 10 A, Aq = 2 meV, Vi = 1 meV, and Umax = 1- We assume Tc = 30 K and 
hence T = 28.24 K. 
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